The Shortest Path Across the Mesoscopic System 
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We study distribution functions (DF) of mesoscopic hopping conductance numerically by search- 
ing for the shortest path. We have found that the distributions obtained by choosing randomly 
the chemical potentials (for a fixed impurity configuration), which corresponds to a typical exper- 
imental situation, coincide with those obtained when both impurity configuration and chemical 
potential is chosen randomly, in agreement with the ergodicity hypothesis. The DFs obtained 
for one-dimensional systems were found to be quite close to the independent predictions of V.I. 
Mel'nikov, A.A.Abrikosov and P. Lee et al. For D — 2, the DFs both for narrow system and thin 
film look similar (and close to the ID case) .The distribution function for the conductance of the 
square sample is nearly Gaussian as predicted by both Altshuler et al and Serota et al. 

PACS numbers: 72.20.Ee, 72.15.Rn, 72.80.Ng, 71.30.+h, 72.15.Eb 



Mesoscopic conductance fluctuations in the insulating 
regime of small, disordered transistors were first observed 
by Pepper in GaAs MESFETs and then studied in de- 
tail in Si MOSFETs by Fowler, Webb and coworkers!) in 
the early 1980s. Extremely strong random fluctuations, 
spanning several orders of magnitude, were observed at 
low temperatures in the conductances of narrow-channel 
devices as the gate voltage was varied. The lognor- 
mal distribution for one-dimensional (ID) conductors 
was long ago analytically predicted by V.I. Mel'nikovjsl 
and A.A.AbrikosovQ. The explanation was provided by 
Lee[^ who proposed a model in which electrons move by 
variable-range hopping (VRH) along a ID chain. Serota, 
Kalia and Lee|^ went on to simulate the ensemble dis- 
tribution of the total chain resistance R and its depen- 
dence on the temperature T and the sample length L. 
In their ensemble, the random impurities are distributed 
uniformly in energy and position along the chain. In ex- 
periments a single device is generally used, so that the 
impurity configuration is fixed, and fluctuations are ob- 
served as a function of some variable external parameter 
such as the chemical potential. An ergodicity hypothesis 
is then invoked to the effect that the same ensemble is 
sampled in both cases, something that has been verified 
experimentally by Orlov et al.^. Using the natural log- 
arithm of the resistance, the authors of Ref. ^ obtained 
for the mean and standard deviation: 
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where ^ is the localization radius and To is the character- 
istic temperature for Mott VRH: Tq ~ l/fc^p^ (p is the 



density of states at the Fermi energy) . It can be seen that 
the size s of the fluctuations decreases extremely slowly 
with length, a result characteristic of ID which was first 
pointed out by Kurkijarvi The explanation is simply 
that exceptionally large resistance elements, even though 
they may be statistically rare, dominate the overall re- 
sistance since they cannot be by-passed in this geometry. 
The averaging assumed in the derivation of Mott's hop- 
ping law for ID does not occur and the total resistance 
takes on the activated form of the largest individual el- 
ement, where ^ is the localization radius and Tq is the 

characteristic temperature for Mott VRH: To = 
(p is the density of states at the Fermi energy) . 

The crossover from ID to 2D were also studied 
pl| . The first theoretical description of hopping conduc- 
tivity in narrow 2D strips was given by Serota|^, and 
numerical simulations from narrow 2D to square 2D at 
a certain tem per ature were also done by X.C. Xie and 
S. Das SarmajlJ. The conductance DF (on the metal- 
lic side) was first fully considered analytically in 2D and 
above by Altshuler, Kravtsov and I. V. Lerner|l^, where 
they have predicted the Gaussian distribution with long 
lognormal tails. Due to the finite widths of channels, the 
2D VRH is also affected by the width of system. To ob- 
serve this effect fully, L. He et al^^ did numerical sim- 
ulations on all 2D cases, including the short 2D which 
has been studied experimentally by Hughes et aZ|]l3j . Al- 
though most of the results are in accordance with these 
of the previous work, the DF for short 2D is of a long tail 
to the low conductance opposed to the low resistance as 
expected. To check the results, we numerically simulate 
these samples by searching for the shortest path across 
system other than by percolating it. The shortest paths 
would be punctures which short out less conductive paths 
in the 2D geometry. The aim of the present work is to 
report this study, and first of all, it is also started by 
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replacing the transport problem across a 2D mesoscopic 
system in Mott hopping regime with a random resistor 
network as did previously. 

Imagine that a particle is about to transport a conduc- 
tive network. In the principle of lowest energy, it prefers 
to the path of the smallest overall resistance among all 
possible paths. The path of lowest energy cost is equiv- 
alent to the shortest path of graph theory, which is just 
what the nature of puncture of the system means. It 
means that we could also approximate the resistance of 
the system by looking for the shortest path, as well as 
by percolating the system as usual. The shortest paths 
among percolation systems have been studied consider- 
ably recently ||l4, 15, 1^, ^ |l^ in terms of minimal path 
or optimal path between a pair of sites within the same 
cluster, and the study on the scaling form of the proba- 
bility of the shortest path with regard to their Euclidean 
distance and the cluster mass (Afa)p^ has shown that 
the average conductance of the percolation backbone is 
strongly correlated with the shortest path, and it de- 
creases with increasing minimal path. This means the 
shortest path determines the average conductance in na- 
ture. 

A network is set up by resistors, py , between sites i 
and j 



Inpy = 2ad + {\E^ - ^| 



\E,-E,\)/2kT, 
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Here a is the inverse localization length, d is the distance 
of two localized sites, Ei and Ej are energies of site 
and fj, is the chemical potential, T is the temperature. 
Energy is chosen randomly from a uniform distribution 
in the range — 0.5 ~ -1-0.5. Thus the mesoscopic system 
is reduced to a random resistor network(RRN). To per- 
colate the network, the resistors joining electrodes are se- 
lected in ascending order until the first percolation path 
connects the reservoirs. The resistance of the percolation 
path is taken to be the resistance of the entire system 

To solve it by the shortest path, Dijkstra algorithm! 19 
is applied, which is efficient in finding shortest paths to 
all nodes from a single source in a fully connected graph. 
In views of graph, the RRN is partially connected, and 
there are more than one source nodes except for the ID 
case. So, some modifications on this algorithm has been 
done in the simulation. 

In calculation, positions of impurities are uniformly 
distributed over the system, their energies are distributed 
evenly between —0.5 ~ -f 0.5, and gate voltage, fi, is ran- 
domly chosen. Thus we can consider the chemical poten- 
tial distributions (for a fixed impurity configuration) and 
the ensemble distributions ( for a fixed chemical poten- 
tial ). For a ID system of L = 1000, chemical potential 
II = and temperature T — 0.001, we have got the pro- 
file of their individual resistances similar to that along 
percolating path(Fig.l of Ref. jl^). It means that it is 
also the single largest hop along the shortest path that 
controls the overall conductance of the system. With 
temperature increasing, the sizes of fluctuation of indi- 
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FIG. 1: The conductance of ID system for v = 0.225: a) 
ensemble distribution function; b) chemical potential distri- 
bution function. The histograms in both figures are the nu- 
merical result 



vidual resistors along the shortest path become close. As 
has been done previously]!^, we also simulate the ID 
systems firstly to validate the simulation. The results 
are shown in Fig.(j^). The size of ID system is 1000 in 
length, and 50 in localization length. The temperature 
is T = 0.001, which gives v = 0.225. The chemical po- 
tential range is /i = —0.1 ~ -1-0.1. As found in Ref.|12|, 
with the increase of temperature, the fluctuation of ID 
resistance become small. 

On the basis of the above validity, the 2D cases are 
numerically simulated by the shortest path, including 
cases of narrow 2D, square 2D and short 2D. In our sim- 
ulation, localization length, temperature, etc., for a 2D 
system, w x L, are the same as in the ID case, while 
w = 100, L = 1000 for narrow 2D, w = 1000 and L = 100 
for short 2D and w = L = 1000 for square 2D system. 
The results are shown in Figj^. From this, it is seen 
clearly that the magnitude of conductance substantially 
increases with the increase of system width, which is in 
agreement with the prediction of both Scrota and Alt- 
shuler et al and the numerical results of X.C. Xie et 
a/pl|] for narrow 2D and square 2D cases. The results 
also reveal that the situations of narrow 2D is close to 
the ID case, and the normal 2D is close to Gaussian. As 
discovered previously [^2], the DF of conductance across 
the short 2D system is still close to that of ID case other 
than a mirror reflection of the latter as expected. Note 
that the peak of short 2D is shifted to the low end of 
resistance more than the other 2D cases. According to 
this, it seems to be more proper to explain the punc- 
ture nature of 2D systems by the position of the peak of 
DF while the long tail of DF shows the bias conductance 
fluctuations. 

In conclusion, the paper has studied the DF of the 
conductance in mesoscopic systems by the shortest path 
which could reflect the puncture nature of VRH system 
more directly. The validity is done and the results are 
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found to be quite close to the numerical results of the 
narrow and square 2D systems by P. A. Lee||, Serota§ 
and X.C. Xie et aZ|l^], i.e., the result for narrow 2D sys- 
tem is close to that oi ID case, and the result for square 
2D is close to Gaussian ||To|. But the unexpected results 
for short 2D system hints that the puncture nature of 2D 
systems could be explained more properly by its peak po- 
sition of DFs and the long tail shows the preference of 
conductance fluctuation. 
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FIG. 2: The ensemble distribution function for the conduc- 
tance of 2D system: a) long system; b) square sample; c) 
short 2D. Histograms are numerical results. The solid curve 
in a) is the RR theoretical fitting for the ID case, and in b) 
is the Gaussian fitting. 
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